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Abstract 

Gauge invariant quark two-point Green's functions defined with path-ordered gluon field 
phase factors along skew-polygonal lines joining the quark to the antiquark are considered. 
Functional relations between Green's functions with different numbers of path segments 
are established. An integral equation is obtained for the Green's function defined with 
a phase factor along a single straight line. The equation implicates an infinite series of 
two-point Green's functions, having an increasing number of path segments; the related 
kernels involve Wilson loops with contours corresponding to the skew-polygonal lines of 
the accompanying Green's function and with functional derivatives along the sides of the 
contours. The series can be viewed as an expansion in terms of the global number of 
the functional derivatives of the Wilson loops. The lowest-order kernel, which involves a 
Wilson loop with two functional derivatives, provides the framework for an approximate 
resolution of the equation. 

PACS numbers: 12.38.Aw, 12.38.Lg. 
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1 Introduction 



Gauge invariant objects are expected to provide a more precise description of observable 
quantities than gauge variant ones. Generally, gauge invariance of multilocal operators is 
ensured with the use of path-ordered phase factors [H [2]. In this respect, the closed loop 
operator, the so-called Wilson loop, showed itself a powerful tool for the investigation 
of the confinement properties of QCD [31 IU [5]. The properties of the Wilson loop were 
studied in detail in a long series of papers P [3 El El HOI HH H2l HH1 HH [15]. 

On the other hand, the usual machinery of quantum field theory, based on the Dyson- 
Schwinger integral equations [161 E] > does not apply in a straightforward way to Green's 
functions of operators involving path-ordered phase factors. The main reason is related to 
the difficulty of obtaining the functional inverses of the nonlocal gauge invariant Green's 
functions and thus of being able to define analogues of proper vertices, which play a 
crucial role in the formulation of integral equations. Expressions of gauge invariant quark- 
antiquark Green's functions in terms of Wilson loops were obtained in the past [HI [191 120] 
with the use of the Feynman-Schwinger representation of the quark propagator [2H [221 (231 
12^] ; these, however, could not be transformed into equivalent integral equations without 
the recourse to approximations related to the quark motion. 

The purpose of the present paper is to investigate the possibilities of deriving integral 
or integro-differential equations for gauge invariant Green's functions which might allow 
for a systematic study of their various properties. We concentrate in this work on the 
quark gauge invariant two-point function, in which the quark and the antiquark fields are 
joined by a path-ordered phase factor, but the methods which we shall develop are readily 
applicable to more general cases. 

Our starting point is a particular representation of the quark propagator in the pres- 
ence of an external gluon field, where it is expressed as a series of terms involving path- 
ordered phase factors along successive straight lines forming generally skew-polygonal 
lines. That representation is a relativistic generalization of the one introduced by Eichten 
and Feinberg in the nonrelativistic case [25]; it was already used in a previous work for 
deriving a bound state equation for quark-antiquark systems [26]; however, in the latter 
work, the bound state equation was derived by circumventing the explicit writing of an 
integral equation for the related Green's function and of the neglected higher-order terms 
of the interaction kernel. One of the main properties of the above representation is that 
in gauge invariant quantities, at each order of the expansion, the paths of the phase fac- 
tors close up to form a Wilson loop. Thus, the corresponding Green's function becomes 
expressed, through a series expansion, in terms, among others, of Wilson loops having 
skew-polygonal contours with an increasing number of sides. 

Several differences occur with respect to the formulation of the Dyson-Schwinger equa- 
tions. First, for the reasons mentioned above, proper vertices are not introduced; instead, 
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we work directly with Green's functions; the various kernels that appear are written ex- 
plicitly in terms of functional derivatives of the logarithm of the Wilson loop average 
and of the quark Green's function. Second, starting from the simplest gauge invariant 
two-point function, constructed with a phase factor along a single straight line joining the 
quark to the antiquark, one generates, through the equations of motion, a chain of new 
gauge invariant two-point functions with phase factors along n-sided skew-polygonal lines 
between the quark and the antiquark (n > 1). On the other hand, every such Green's 
function (with n skew-polygonal sides) can be related with the aid of functional relations 
to the lowest-order Green's function (n = 1) and thus, in principle, an equation involving 
only the latter Green's function is possible to construct. The third difference arises at 
the level of the presence of nested kernels, which do not occur in the Dyson-Schwinger 
equations and which persist here due to background effects induced by the Wilson loops: 
each nested kernel is modified by its new background when inserted inside a higher-order 
Wilson loop as compared to its original expression. The remaining terms in the kernels 
have the property of conventional irreducibility. 

The integral equation that we obtain is constructed as an expansion in terms of the 
global number of derivatives of the logarithm of the Wilson loop average. Although it 
involves an infinite series of kernels and Green's functions, at each order of the expansion 
the explicit expressions of the kernels and of the relations between high-order Green's 
functions with the lowest-order one can be obtained from definite formulas. 

On practical grounds, an increasing number of derivatives of a Wilson loop, each 
derivative occurring on a different region of the contour, is generally expected to give a 
relatively decreasing contribution at short- and at large-distances. Therefore, the series 
expansion of the kernels in terms of functional derivatives of Wilson loops can also be 
considered as a perturbative expansion, the most important contribution coming from 
the lowest-order non-vanishing term and involving the smallest number of derivatives. 
That property allows us to consider solving the integral equation with appropriate ap- 
proximations. 

The plan of the paper is the following. In Sec. 2, we introduce the definitions and 
conventions that will be used throughout this work. Section 3 deals with the representa- 
tion of the quark propagator in external field in terms of path-ordered phase factors. In 
Sec. 4, functional relations are established between various Green's functions. In Sec. 5, 
the integral equation for the gauge invariant quark two-point function with a straight line 
path is established and the structure of the kernel terms is displayed. Section 6 deals with 
the question of analyticity properties of the Green's function. A summary and comments 
follow in Sec 7. Two appendices are devoted to the presentation of the summation method 
with free propagators and the study of the self-energy function. 
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2 Definitions and conventions 



We introduce in this section the main definitions and conventions that we shall use 
throughout this work. 

We consider a path-ordered phase factor along a line C yx joining a point a; to a point 
y, with an orientation defined from x to y: 

-ig [ V dz^AJz) 

U(C yx ;y,x) = U(y,x) = Pe h , (2.1) 

where = J2 a A a J,° A a ^ (a = N% — 1) being the gluon fields and t a the generators 

of the gauge group SU(N C ) in the fundamental representation, with the normalization 
trt a t b = \$ ab ■ A more detailed definition of U is given by the series expansion in the 
coupling constant g\ all equations involving U can be obtained from the latter expression. 
Parametrizing the line C with a parameter A, C = {x(A)}, < A < 1, such that x(0) = x 
and x(l) — y, a variation of C induces the following variation of U [U(x(X),x(X')) = 
U(X,X%A{x(X))=A(X)): 

#7(1,0) = -ig5x a (l)A a (l)U(l,0) + igU (1,0) A a (0) 5 x a (0) 

+ig f 1 dXU(l, X)x ,f3 (X)Fp a (X)5x a (X)U(X,0), (2.2) 



o 

where x' — || and F is the field strength, F^ v = d iii A v — d u A^ + ig[A^, A v \ . The variations 
inside the integral lead to functional differentiation of U, while the variations at the end 
points (marked points) are defined as leading to ordinary differentiation. The functional 
derivative of U with respect to x(A) (0 < A < 1) is then [6]: 

= igU(l, A) x' p (X)Ff} a (A) U (A, 0). (2.3) 

For paths defined along rigid lines, the variations inside the integral in Eq. (12.21) are 
related, with appropriate weight factors, to those of the end points; a displacement of one 
end point generates a displacement of the whole line with the other end point left fixed. 
Considering now a rigid straight line between x and y, an ordinary derivation at the end 
points yields: 

9U Qy a X) = -igMvMy, x) + %g(y - xf £ dX A U(l, X)Fp a (X)U(X, 0), (2.4) 

dU ^ a X) = +igU{y, x)A a (x) + ig(y - xf £ dX (1 - A) U(l, X)Fp a (X)U(\, 0). 

(2.5) 

When considering path variations of gauge invariant quantities, with paths made of 
segments, the end point contributions involving the explicit A a terms disappear, being 
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cancelled by similar contributions coming from neighboring segments or from variations 
of neighboring fields. The general contributions that remain at the end are those coming 
from the internal part of the segments represented by the integrals in Eqs. (I2.4p - (I2.5D . 
We adopt the following conventions to represent such contributions: 

5 -B^L = ig(y - xf £ dX A U(l, X)Fp a (X)U(\, 0), (2.6) 

- ig(y - xf f 1 d\ (1 - A) 17(1, X)F /3a (X)U(X, 0). (2.7) 



5x a - Jo 

The first equation above corresponds to a displacement of the end point of the segment 
(taking into account the orientation on the path), while the second equation corresponds 
to a displacement of the starting point of the segment. Equations (12.41) and (12.51) can be 
written as 

«^ = ^ W(y ,, ) + ^, (2.8) 

^ = +igU(vM) + ^. (2.9) 

If two phase factors U(z, y) and U(y, x) along segments are joined at the point y (a 
marked point), then, with the aid of the previous notations, we have: 

S ( TII >„, >\ SU(z,y) Ir , , , „, JU(y,x) 



^ {Uiz , y)U{y , x) ) = ^ U{y , x) + U{z ,yf^>. ( 2 ,0) 

The Wilson loop, denoted 3>(C), is defined as the trace in color space of the path- 
ordered phase factor ( 12. lft along a closed contour C: 

1 -iq <p dx^AJx) 
$(C) = — trPe "Jo \ (2.11) 

where the factor 1/N C has been put for normalization. It is a gauge invariant quantity. 
Its vacuum expectation value is denoted W(C): 

W(C) = ($(C)>, (2.12) 

the averaging being defined in the path integral formalism. 

We shall represent the Wilson loop average as an exponential function, whose argument 
is a functional of the contour C [H [14"] : 

W{C) = e F{ ~ C \ (2.13) 

In perturbation theory, F(C) is given by the sum of all connected diagrams, the connection 
being defined with respect to the contour C, after subtraction of reducible parts [T3] , 
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Variations of W(C) due to local deformations of the contour C can then be expressed in 
terms of variations of F(C): 



SW(C) 



5x c 



SF(C) 



W(C). (2.14) 



This property is also generalized to the case of rigid variations of paths (segments). If the 
contour C is a skew-polygon C n with n sides and n successive marked points xi, x 2 , . . ., 
x n at the cusps, then we write: 

W{x n , x n _i, ..., Xl ) = W n = e Fn ( Xn > • • • ' Xl ) = e Fn , (2.15) 

the orientation of the contour going from x\ to x n through x%, X3, etc. (i.e., towards 
xs with indices increasing by one unit). Then, according to the definitions fl2.6l) - fl2.10l) . 
the notation 8F n /8x^ means that the derivation acts on the internal part of the segment 
XiXi + \ with Xi + i held fixed (x n+ i = x±), while SF n /5xf means that the derivation acts on 
the internal part of the segment Xi-iXi with a?j_i held fixed (xo = x n ). 
The gauge invariant two-point quark Green's function is defined as 

S aj3 (x,x']C x > x ) = — — - {^p{x') U(C x > x ;x',x) i^ a {x)), (2.16) 

a and (3 being the Dirac spinor indices, while the color indices are implicitly summed. In 
the present work we shall mainly deal with paths along skew-polygonal lines. For such 
lines with n sides and n — 1 junction points yi, 1/2, ■ ■ ., y n -\ between the segments, we 
define: 

S( n )(x,x';y n -.x, ..., yi ) = -— (if>(x')U(x' , y n -i)U(y n - 1} y n - 2 ) ■ ■ ■ U(y u x)if>(x)). (2.17) 

The simplest such function corresponds to n — 1, for which the points x and x' are joined 
by a single straight line: 

S {1) (x,x') = S(x,x') = -— ($(x')U{x',x)^{x)). (2.18) 

(We shall generally omit the index 1 from that function.) 

The free propagator will be designated by So (without color group content): 

S (x,x') = S (x-x')= /-^P e-W*-*') (2.19) 

J (27r) 4 7-p — m + is 

In conjunction with the definitions fl2.6p - fl2.7l) . we shall also introduce the notations 

fci) = -— ®{x>)~ 5U - {x '' x) v>(rr)), Zs{x ; x ' ] = -- (W)~ 5U - {x r x) ^)). 

Sx»- N c WK ' 5xv- VK }h 5x'»+ N c WK ' 5x'» + VK " 

(2.20) 
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3 The quark propagator in external field 



We shall use a two-step quantization method, by first integrating the quark fields and 
then, at a second stage, integrating the gluon fields through Wilson loops. The first 
operation yields among various quantities the quark propagator in the presence of an 
arbitrary external gluon field. The latter, designated by S(x,x'; A), or by S(A) for short, 
satisfies the usual equation 

(i'y.d^ — m — g"Y.A(x)^jS(x, x'; A) = i5 4 (x — x'). (3.1) 

To exhibit a Wilson loop structure in gauge invariant quantities, it is necessary to describe 
the quark propagator in external field by means of path-ordered phase factors. To this 
end, we shall first introduce a representation, already used in Ref. [26], which combines 
path-ordered phase factors along straight lines and free quark propagators. At a later 
stage, to sum directly self-energy effects, we shall replace the free quark propagator by 
the full gauge invariant Green's fuction (12.181) . 

The starting point of the representation is the gauge covariant composite object, de- 
noted Sq(x, x'), made of a free fermion propagator So(x, x') (without color group content) 
multiplied by the path-ordered phase factor U(x, x') [Eq. (12. ip ] taken along the straight 
line xx 1 : 

'S (x,x')] a b = S (x,x')[u(x,x')] a b . (3.2) 

[a, b: color indices.] The advantage of the straight line over other types of line is that 
under Lorentz transformations it remains form invariant and in the limit x' — > x U tends 
to unity in an unambiguous way. S satisfies the following equation with respect to x: 

(i r d {x) -m- gj.A(x)) S (x, x') = i5\x - x') + i ^j^ S {x, x'). (3.3) 

A similar equation also holds with respect to x', with x held fixed, with the Dirac and 
color group matrices acting from the right. 

The quantity — i(ij.d/ x \ —m — g'y.A(x))S A (x — x') is the inverse of the quark propagator 
S(x, x'; A) in the presence of the external gluon field A. Reversing Eq. (13. 3ft with respect 
to S(A)~ 1 , one obtains an equation for S(A) in terms of Sq: 

S(x, x'; A) = S (x, x 1 ) - J d'x" S(x, x"; A) 7 Q ^ff^ ■ (3.4) 

Using the equation with x', or making in Eq. (13 .4p an integration by parts, one obtains 
another equivalent equation: 

S{x, x'- A) = S {x, x') + J d 4 x" 8S f^*P T S{x", x'- A). (3.5) 

Equations (13.41) or (13. 5p allow us to obtain the propagator S(A) as an iteration series 
with respect to So, which contains the free fermion propagator, by maintaining at each 
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order of the iteration its gauge covariance property. For instance, the expansion of Eq. 
(13.41) takes the form: 



S{x,x';A) = S {x,x') - / d 4 yiS'o(x,yi)7 



+ / d 4 yid 4 y 2 S (x,y 1 )~f' 



n SS (y 1 ,y 2 ) a2 5So(y 2 ,x' 
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+ ■ ■ 



(3.6) 



5yT vy* 

Equations (13 .4p and (13. 5p are relativistic generalizations of the representation used for 
heavy quark propagators starting from the static case [2"B] . 

In order to sum, for later purposes, self-energy effects, one can use for the expansion 
of the propagator S(A), instead of the free propagator S , the full gauge invariant Green's 
function (12. 18p . To this end, we define a generalized version of the gauge covariant object 
S [Eq. <ET2} ]. by replacing in it S with S [Eq. (EH])]: 



S(x,x') b = S(x,x') U(x,x f ) 

The Green's function S satisfies the following equations of motion: 

5S(x, x') 



(i~f.d( x ) — m)S(x, x') = i5 (x — x') + iY 



5x^~ 



S(x,x')H 7 - d&) -m)=i5\x-x')-i ^ S j^ ) Y- 



(3.7) 

(3.8) 
(3.9) 



[Notice that the orientation of the path in S(x, x') is from x to x'.] Then S satisfies the 
equation 



(Vy.<9( x .) — m — g'y.A(x)^ S(x, x') = id 4 (x — x) 



a /SS(x,x') 



U(x,x') + S(x,x') 



, 5U(x,x'] 



5x a ~ 8x a + 
from which one deduces the expansion of S(A) around S: 



(3.10) 



S(x, x'\ A) = S(x, x')U(x, x') - S(x, y; A) 7 C 



8S(y,x') 5U(y,x') 
— = U(y,x ) + b(y,x ) — = — ; — 



(3.11) 



[The integrations on intermediate points are implicit.] Using the equations of 5* and 
S(A) relative to x', or making in Eq. (13. lift an integration by parts, one obtains another 
equivalent equation: 



S(x,x'; A) = S(x,x')U(x,x') + 



SS{x,y) 



U(x,y) + S(x,y) 



SU(x,y) 



YSiy^x^A). (3.12) 



5y a+ 5y a ~ 

A graphical representation of Eq. (13.121) is shown in Fig. [1] 

Equations (13. 1H) - (13.1 2D constitute the basic formulas that will be used to express 
equations of motion of gauge invariant quark Green's functions in terms of Wilson loops 
and gauge invariant two-point Green's functions. 



S 



y v 




Figure 1: Graphical representation of Eq. f 1 3 . 1 2 j) . The double line (one full and one dotted 
joining two circles) represents the gauge invariant Green's function S(\) = S [Eq. (12.181) ] 
with a path along a single straight line; the single dotted line represents the phase factor, 
the dashed line the quark propagator in the external gluon field, the arrow the orientation 
on the path. The cross represents the rigid path derivation with one end fixed; it is placed 
near the end point which is submitted to derivation, y is an integration variable. 



4 Functional relations for Green's functions 

Functional relations between various gauge invariant quark Green's functions are obtained 
with a systematic use of Eqs. (13. lip or (13.121) . 

Let us consider the Green's function S( n ) [Eq. ( 12. 17ft ]. Integrating with respect to the 
quark fields, one obtains: 

S {n) (x, x'; y n _i, . . . , y x ) = -^r (U (x', y n -i)U(y n - l7 y n - 2 ) ■■•U (j/i, x)S(x, x'\ A)). (4.1) 

The simplest case of this equation, corresponding to n — 1, is: 

S m {x,x') = S(x,x') = ±- (U(x',x)S(x,x'; A)). (4.2) 

The quark field integration yields also a corresponding determinant, which is a functional 
of the quark propagator in the external gluon field A. That determinant will not, however, 
play an active role in the subsequent calculations and hence will not explicitly appear in 
the various formulas that we shall meet; it will rather contribute as a background effect; 
in particular, it contributes to the evaluation of the Wilson loop averages, unless the 
quenched approximation is adopted. Therefore, the averaging formulas that we shall 
encounter should be understood with the presence of the quark field determinant. The 
expansions that will be used for the quark propagator in the external gluon field can also 
be repeated inside the quark field determinant if Wilson loop averages are to be evaluated. 
Using now for S(A) Eq. (I3.12p . one obtains: 

S( n ) (x, x'\ y n _i, . . . , y x ) = j^S(x, x') (U (x', j/ n _i) ■■■U(y 1 , x)U(x, x')) 



■±(U{a/, y n -i) ■ ■ ■ Ufa, x) {^jh^U(x, y n ) + S(x, y n ) 5 -^^ 
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S(x, x') e Fn +^ x '' yn ~ 1 ' ■■■yy^ x ) 



+ 



5S(x,y Tl 



5 



+ S ( x 'Vn)7-^) ^° S (n+i)(yn,x';y n -u ...,y h x). 



A graphical representation of this equation for n = 3 is shown in Fig. [2j 



(4.3) 



Hi 



yi 



Hi 



J/2 



Hi 



1)2 



11 4 



111 



>(3) 



3/3 S, 



. il l 



(4) 



Figure 2: Graphical representation of Eq. (14. 3 p for n = 3. Same conventions as in Fig. [U 
?/3 is an integration variable. 



The Green's function Si n \ satisfies the following equation of motion with x: 



(ij-d( x ) - m)S {n) (x,x';y n _i, . . . ,y x ) = i5 A (x - x')e 



F n (x,y n _i,...,yi) 



+«7 



p&S^fog'iS/n-i, . . . ,yx) 



(4.4) 



A graphical representation of this equation for n — 1 and n = 3 is shown in Fig. [3J 



(17. 9* - m) 9- 



a; S 



id 4 (x - x') 



(i-y.d x - m) . • = V + ^ # 

x 5 C 3) ij 4 (a;-a/) x S{3) x > 

Figure 3: Graphical representation of the equations of motion of S'(i) = S and S@). Same 
conventions as in Fig. dJ 
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5 Integral equation 



The equations of motion of the gauge invariant Green's functions S(n) [Eqs. (13. 8ft and (14.411 
involve in their right-hand sides as unknowns the rigid path derivative of the Green's 
functions. The core of the problem amounts therefore to the evaluation of the rigid 
path derivative of Green's functions. That task, however, is facilitated by the functional 
relations (14. 3p . which relate two successive Green's functions with increasing index. They 
allow the evaluation of the rigid path derivative of a Green's function in terms of a similar 
derivative of a Wilson loop average and the derivative of a Green's function with a higher 
index. Systematic repetition of this procedure allows us therefore to express the rigid path 
derivative of a Green's function in terms of a series of Green's functions whose coefficients 
are functional derivatives of Wilson loop averages. One thus obtains chains of coupled 
integral (or integro-differential) equations between the various Green's functions. At the 
end, each Green's function can be expressed, at leading order of an expansion, by 
means of the functional relation (14. 3p . in terms of the lowest-order Green's function S, and 
thus an equation where solely the Green's function S would appear becomes reachable. 

In the present work we are mainly interested by the simplest Green's function S 
and therefore we shall concentrate our considerations on the equation of motion of that 
quantity. 

The rigid path derivative of Sr n ) along the segment xy\ is obtained from Eq. (]4.3j) : 

^(x,^;^,...,;/!) = e F n+1 (x', y n . u ..., Vl ,x) g, n 

5xf- 5x»- ' 

+ { Ey „+ + S(x, Vn)j^) 7 S (n+1) (y n , x ; y^, ...,y x ,x). (5.1) 

Eliminating in the right-hand side of the latter equation the product e Fn+1 S through Eq. 
(14.31) . one obtains the equation 

8S { n){x,x'\y n _ u ...,y 1 ) _ 5F n+1 (x', . . . , y u x) 



5x^~ Sx 11 ' 

/ 6 SF n+1 ^ ( 5S{x,y n ) 
+ (^r " j^r) {-^ y¥ - + S(x, y n )—) 7 S {n+1) (y n , x ; y n . u ...,y u x). 

(5.2) 

[Integrations on new variables in the right-hand sides are implicit.] For n — 1, one has: 

5S(x,x') SF 2 {x',x) , 
— r = — ? S(x,x) 

We next evaluate, in Eq. (15. 3p . the action of the path derivation operators on 5(2). 
The operator 8/8x^~ acts here on the segment xx' and therefore it can be brought without 



S( n )(x,x';y n -i, ...,yi) 
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harm to the utmost right, where an equation similar to Eq. (15. 2p (with a relabelling of 
some arguments, the point x being now a junction point on the path y^xx' of S7 2 )) is used 
with n = 2 and then 8F 3 /8x^~ is brought back to the left; during the last operation it 
is also submitted to the action of the operator 8/8y ai ~ . The resulting terms that involve 
5(2) are: 

8F 3 (x',x, yi ) 8F 2 (x',x)^ (5S(x :yi ) 8 s , 

+ 6 7x^ Vl) Vl) 7 " S{2) x ' x) - (5 ' 4) 

Next, one observes that 8S^/8y ai ~ is part of the equation of motion of 5(2) [Eq. (|4.4[) ] . 
yx being now one of the fermionic ends of Sm- Using the latter equation and making an 
integration by parts with respect to yi, one arrives at a simplified expression in which 
in the second derivative of F 3 the derivation S/Sy" 1- , which is along the segment yix, 
is replaced by the derivation 8/8yi 1+ , which is along y\x' . [The delta functions of the 
equations of motion do not contribute here because of the existence of the difference term 
{8F 3 /8x»~ -8F 2 /8x^).\ 

The net result is, including also the resulting 5( 3 ) terms: 

8S(x,x') 8F 2 {x',x) , 8 2 F 3 {x',x,y 1 ) 

= S{x, x ) - — -p S(x, Vl ) 7 S {2) ( yi , x ; x) 

>5S(x,yi) t 8 \ ( 8 8F 2 (x,x') 



x ( ^f } + S(y u V2)j^=) 7° 2 S {3) (y 2 , x'; x, Vl ). (5.5) 

In the term containing 5(3) , the factors with the derivatives with respect to x and 
y 2 are treated in the same way as were those with x and y\ with Sr 2 y, they yield at the 
end the factor 8 2 F 4 /8x^~8y 2 2+ plus a term with 5(4) having a similar structure than the 
one with S7 3 ) above. Repeated use of the procedure described with 5(2), yields a series 
expansion in 5( n ) (n = 2, 3, . . .) where all terms have similar structures. One obtains: 

5S(x,x') 5F 2 (x',x) , 8 2 F 3 (x',x iyi ) 

= S(*, - ) ~ - Sx ,_- W+ S(x, y x ) 7 S {2) (y u x ; x) 

(5S(x, yi ) 8 \ 8 2 F 4 (x' J x J y u y 2 ) , 

-{ x 7;Ql+ +S{x, yi ) J —) 1 . . + S{y u y 2 ) 1 >S {3) (y 2 ,x;x, yi ) 



8y^ + v 7a± '6y?-' ' 8x^8y° 2 



> 8S[x, yi ] 
5y? 1+ ' ~ v ~ 
/8S(y n _ 3 ,y n ^ 2 ) 

w 8 -F 1 n +i(x , x, yi, . . . , y n -i) , \ a „-i or > \ in r\ 

x = - an _ 1+ 5(y n _2,Z/n-i)7 1 5 (n) (y n _i, x ;x, yi,..., y n _ 2 ). (5.6) 



Ctn-2 
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Equation (15. 6p determines the action of the rigid path derivative on the Green's func- 
tion. It is expressed in terms of derivatives of logarithms of Wilson loop averages appearing 
in a series of terms with skew-polygonal type contours. Although the resulting expression 
still contains other derivatives along internal lines, as well as implicit terms, their action 
will be determined with respect to the structure defined by Eq. (15.61) . The derivatives of 
the logarithms of Wilson loop averages, together with the accompanying quark Green's 
function S, play here the role of kernels of the integral equation we are searching for; they 
are the analogs of the kernels made of propagators in Feynman diagrams appearing in 
Dyson-Schwinger equations. By analogy, we shall often call them diagrams. 

We now study the structure of the kernels that are present in the expansion (15.61) . 
We notice that in the term with S( n ), the utmost left derivative related to x is connected 
to the utmost right derivative related to y n -\ through the term 6 2 F n+ i/6x~Sy^_i] this 
does not leave room for the existence of reducible type terms made of disjoint subsets of 
connnections; such terms are actually parts of the definitions of the S( n )S when expanded 
in terms of free propagators. In the present case, all remaining derivatives 5/5yi (i = 
1, . . . ,y n -2) either will act within the skew-polygonal line xy\ . . .y n -\ or will leave that 
line to be connected to a larger contour associated with an S( m ) with m > n. 

To exhibit the latter feature, we consider in Eq. ( 15.61) the general term St n ) and its as- 
sociated derivative terms. We notice that the derivatives S/5yi {i = 1, . . . , y n -2) no longer 
act on the fermionic end of Sr n ) (represented by y n -i) but only on the internal junction 
points of the skew-polygonal line xy\ . . . y n -i, hence, one cannot use the equation of motion 
of S( n ) with these variables. We first consider the action of S/5y~_ 2 . It is brought to the 
right to act on Sr n y, during this operation it may also act on the term 8 2 F n+ i/8x~8y^_ 1 
to yield a third-order derivative of F n+ i. When in front of S( n ), it is convenient, to 
maintain symmetry with S/Ey^-i, to replace 8/5y~_ 2 by (<9/ch/„_ 2 — 8/5y^_ 2 ); the total 
derivative is then used for an integration by parts to convert 5 3 F n+ i/5x~5y^i5yn-2 into 
S 3 F n+ i/ 'Sx~Sy^_ 1 Sy^_ 2 - The derivative S/5y^- 2 acting on S( n ) has two kinds of effect, 
according to Eq. (14. 31) (with a relabelling of variables): in the first place, it yields the 
derivative term 5F n+ i/5y^- 2 and in the second, it generates a new Green's function S( n +i) 
with the fermionic end point y n (the global result is very similar to that of Eq. (15.21) with 
a relabelling of variables). The term 5F n+ \/ '5y^_ 2 is an insertion along the line xy\ . . . y n -i 
and may also be submitted to other derivations coming from the remaining variables yj 
(j — 1, . . . , n — 3). In the newly generated term with Sr n +u, the actions of the deriva- 
tives 5 /8y^_ 2 and S/5y~ can be combined to give 5 2 F n+2 /5yn- 2 Sy^ . This term represents 
now a connection between the line xy\ . . . y n -i and the segment y n -\yn\ it crosses the 
connection line xy n -\ which had already appeared with the term 5 2 F n+ x/5x~5y^_ 1 ; the 
two connections can therefore be interpreted as forming a crossed diagram. We observe 
that the appearance of a crossed diagram has been accompanied with the increase of the 
number of segments of the contour by one unit. 
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The action of the remaining derivatives S/SyJ (j = 1, . . . , n — 3) can be studied in a 
similar way. The qualitative features displayed up to now remain unchanged. 

We can summarize the above results by grouping the terms that appear in front of a 
Green's function S( n ) into three categories: (i) A term that is completely connected: F n+ i 
is submitted to the n derivations 5 n /5x~5yf . . . 5y^-i- (h) Crossed diagrams that involve 
at least one F n+ i and some other Fs with lower indices, (iii) Nested diagrams, represented 
by insertions within the connection line xyi . . . y n -i or within smaller connections of that 
line or within crossed diagrams. As a general property, no terms of reducible type (disjoint 
connections) exist. 

The general structure of the derivative 5S/5x^~ is: 

S(x,x') + K 2 ^-(x',x,y 1 ) S {2) (y 1 , x'\ x) 

oo 

+ J2K ilM _(x',x,y 1 , . . . , 2/i-i) S^iyi-i^x'-x^y!, . . . ,yi- 2 ), (5.7) 

i=3 

where the kernels Ki (i — 1, . . . , oo) are composed of the three categories of terms quoted 
above and of (i — 1) quark propagators S, and eventually of their derivatives, along the 
segments of the {i + l)-sided skew-polygons. The total number of derivatives contained 
in K n is n. 

The explicit expression of 5S/5x^~ up to the fourth-order of its expansion is: 



5S(x, x') 
5x»~ 



SS(x,x') 8F 2 (x',x) a . 8 2 F 3 (x',x, yi ) 



5x^~ Sx 11- 
<5 3 F 4 (x',£, y u y 2 ) 



S(x, x') 



Sx^Sy? 1 - 



S{x,yi)"i ai S^{y x ,x']x) 



+ 



S(x, y x ) 7 Ql S(y u y 2 ) 7 02 S {3) (y 2 , x'; x, y x ) 



5x^-5y^ 1+ 5y^ + 

4 S(x, Vl ) 7 - + S(y u y 2 )^) 7 - S (3) (y 2 , x>- x, Vl ) 



5 2 F, 



5x^~5y 2 



5*F 5 (x', x, yi, y 2 , y 3 ) 5 2 F 4 (x', x, y u y 2 ) 5 2 F 5 (x', x, y u y 2 , y 3 ) 



5x^Sy^ 1+ Sy^ + Sy^ + Sx"-Sy^ 2+ 
x S(x, yx) 7 Ql S(y u y 2 ) ^ S{y 2 , y 3 ) 7° 3 S {4) (y 3 , x'\ x, y u y 2 



S{x, Vl )^ 
S 3 F 5 



6 3 F, 



5x»-5yr + 5y c 3 



5S(y u y 2 ) 



S{yi, y%) 



SF 5 
SF 5 



5x»-5y^5yr + S(Vl > ^ 7 " + S ^ ^ 5y? 

f S S(y 2 ,y 3 ) 



5 2 F, 



6 F* 



5xn-5y^ + V 5y{ 



5 2 F, 



Sx^Sy^ + Sy^ 1+ Sy^ 



fSS(y 1 ,y 2 ) 

S(yi,y2,)l a2 S(y 2 ,y 3 ) 



7 a2 S(y 2 ,y 3 ) 



S{y 2 ,y 3 )Y7^+ 



'2 Sy 2 
l a ' A S(4)(y3,x';x,y u y 2 ) 



(5. 
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The expansion, up to third-order terms (including 5(3)), is represented graphically in Fig. 

m 




Figure 4: The expansion of SS/5x up to third-order terms. Same conventions as in Fig. 
[TJ The ys are integration variables. 

In Eq. (I5.8p . the term in front of Sp) and the first terms in front of 5( 3 ) and 5(4) 
correspond to the completely connected diagrams. The second term (within the first 
brackets) in front of Su\ represents a crossed diagram. The remaining terms with 5( 3 ) 
and 5(4) correspond to nested diagrams. 

The persistence in the integral equation of nested diagrams may seem puzzling. In the 
Dyson-Schwinger equation for the self-energy, once the internal propagators of diagrams 
are replaced by full propagators, no nested diagram survives; this is equivalent to stating 
that the summation of all internal diagrams in nested diagrams yields the full propagators. 
In the present formalism, the survival of nested diagrams is a consequence of a background 
effect induced by the Wilson loop; each sub-diagram is actually calculated in the presence 
of the Wilson loop appearing at the order of the global diagram. At each order of the 
expansion the contour of the relevant Wilson loop changes with an increase of the number 
of segments forming the contour. The insertion of a low-order diagram in a higher-order 
diagram thus modifies the expression of the former, since it is now expressed with the 
new Wilson loop derivatives. In dealing with ordinary Feynman diagrams, one does not 
encounter the above background generated effects. 

These statements can be explicitly checked in Eq. (15.81) . A way of ignoring background 
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effects can proceed as follows: (i) Assimilate an n th -order derivative of a function F m 
(m > n+1) to an ordianary n-point function with propagators attached to the end points 
of the corresponding segments, with the convention that a derivation of the type 5/Sy~ 
can be converted into a derivation of the type 5 / 5y + with a change of sign, (ii) Consider 
first the approximation where S( n ) is given by the first term of the expansion (14. 3j) . (iii) 
When two loop contours corresponding to Wi and Wj have a common segment, replace 
the whole by a single contour corresponding now to Wi+j-2. It can then be verified that 
a systematic expansion of the terms 5S/5y~ inside the nested diagrams with an iterative 
procedure cancels order-by-order all other terms present in the nested diagrams and thus 
the latter disappear, as expected, when no background effects are retained. At a second 
stage, considering the further terms of the expression of S( n ) [Eq. (14. 3p ] one finds that 
the latter are themselves background generated effects and the repetition of the above 
procedure makes them in turn disappear. 

The rules of appearance of nested diagrams can be deduced from Eq. (15.61) and verified 
in Eq. fl5"B . 

In two dimensions, where the Wilson loop averages are determined by the areas of the 
surfaces lying inside the contours [13], the second functional derivatives of the functions 
F are delta-functions and in general the nested diagrams disappear. One might expect 
that, in four dimensions, the nested diagrams, even if not disappearing, remain negligible 
on quantitative grounds. A more complete idea of their role could be obtained only when 
renormalization properties of the integral equation are studied. 

Equation ( 15. 6ft . together with relations ( 15.21) . allows the calculation of the term SS/Sx~ 
through an expansion involving an increasing number of functional derivatives of Wilson 
loops. The calculation of the expression of the kernel K n appearing in the expansion (15. 7p 
requires solely consideration of terms of order lower or equal to n. The integral form of 
the equation of motion ( 13.81) is: 

S(x,x') = S (x,x') + J d 4 x"S (x,x")j» 6S j*''^'\ (5.9) 

in which one has to inject the expression of 5S/5x~ resulting from Eq. (15. 71) . 

At short- distances, governed by perturbation theory, each derivation introduces a new 
power of the coupling constant and therefore the dominant terms in the expansion are 
the lowest-order ones. At large-distances, Wilson loops are saturated by the minimal 
surfaces having as supports the contours El [26]. Here also, the dominant contributions 
come from the lowest-order derivative terms. Therefore the expansion in Eq. ( 15.71) can be 
considered in general as a perturbative one whatever the distances are, provided that for 
each type of region the appropriate expressions are used for the Wilson loops. The first 
term of the expansion, represented by a single derivative, is null for symmetry reasons (the 
derivative 5F2(x', x)/6x~ being orthogonal to xx'). Hence the non-zero dominant term 
of the expansion is the second-order derivative term. Furthermore, the various Green's 
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functions S( n \ are themselves dominated by their lowest-order expression of Eq. (14.31) . 
involving only S and a Wilson loop. In that approximation, 5S(x, x')/5x~ takes the form 



Thus, the dominant part of the integral equation (I5.9P to be solved reduces to a closed 
form involving only the full propagator S, the free propagator So, a Wilson loop average 
and its second-order rigid path derivative. 

The lack of manifest symmetry in Eq. (15.91) between the coordinates x and x' is due 
to the presence of the closed contours of the Wilson loops, which do not allow immediate 
factorization of propagators through convolution operations. Nevertheless, because of 
translation invariance, the Green's function S(x, x') depends only on the difference (x—x f ); 
therefore, once the integrations in Eq. ( 15. 9ft are done, one should recover the desired 
symmetry. If, instead of the equation of motion (13.81) . relative to x, we had used the 
equation of motion (13.91) . relative to x' , we would have found an integral equation where 
So(x",x') acts on 5S(x, x")/8x from the right. 

As a complementary remark with respect to the method of approach developed in 
this work, we point out that another way of proceeding would consist of trying to sum 
diagrams constructed with free quark propagators associated with phase factors. This 
method actually corresponds to the use of the expansion equations (13.41) and (13.51) . instead 
of ( 13. lip and ( 13.121) . It, however, becomes rapidly intricate due to the presence of the 
background Wilson loop effects. Nevertheless, the first few terms of Eqs. (15. 7ft and ( 15. 8ft 
can be reconstituted rather easily. The main aspects of this method are presented in 
Appendix [A] 

The integral equation ( 15.91) . together with the expressions (15. 7ft or (15.81) . did not make 
any explicit reference to the quark self-energy function. The latter can be obtained once 
the Green's function S is calculated, through its inverse. It is, however, also possible to 
construct it directly, by setting up a specific equation for it. This is presented in Appendix 



The problem of solving Eq. (13.81) . or its equivalent Eq. ( 15.91) . together with expression 
( 15.71) . leads us to the question of representation of the gauge invariant two-point Green's 
function. This will be considered in Sec. |6j 

6 Analyticity properties of the Green's function 

One of the advantages of paths along straight lines is the fact that the expressions of the 
corresponding Green's functions become dependent only on the end points of the paths. 
This feature in turn allows a simple transition to momentum space by Fourier transfor- 
mation. Much of the informations on Green's functions are provided from momentum 
space, since it is there that their spectral properties are determined. 



8S(x, x') 



5xv- 




e F 3 (x', x, yi) s ^ y j 7 «! 3 ( yij x ,y ( 5 1Q ) 



m 
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From this point of view, the quark two-point gauge invariant Green's functions hold 
a particular position. Because of confinement of colored objects, it is not possible to cut 
the path joining the quark to the antiquark by inserting in it a complete set of physical 
states, which are color singlets. This feature seems to suggest that gauge invariant two- 
point Green's functions should not have any singularities. 

The situation is, however, more complex than it seems. Gauge invariant two-point 
Green's functions possess singularities originated from perturbation theory. This is cor- 
roborated by the integral equation (15.91) . in which the presence of the free quark prop- 
agator generates new singularities in the complete solution. An analysis, starting from 
perturbation theory, is therefore necessary. 

We shall admit that, in a domain where perturbation theory is valid, it is meaningful 
to consider quarks and gluons as physical particles with positive energies, described by 
corresponding physical states. It is then advantageous to consider the path-ordered phase 
factor (12. ip in its representation given by the series expansion in terms of the coupling 
constant g, the n th -order term of the expansion containing (n — 1) gluon fields (n > 1). 

Adopting here an operator formalism, we observe that the gauge invariant quark two- 
point function involves two kinds of orderings for its defining fields. The first is the path- 
ordering (or P-ordering) which concerns the color index arrangements of the gluon fields 
according to their positions on the path. The second is the time-ordering (T-ordering) 
or chronological product which enters in the definitions of Green's functions and operates 
once the P-ordering is done. 

Another advantage of paths along straight lines is that once the timelike or spacelike 
nature of the distance between the quark and the antiquark is fixed, the nature of the 
mutual distances of the gluon fields in the Green's function S [Eq. (I2.18p ] is also fixed 
in the same way, because of their alignement along the segment joining the quark to the 
antiquark. Therefore, the chronological product of the n th -order terms in S reduces to 
two terms, defined by the relative time between the quark and the antiquark. According 
to the definitions (12. ip and (I2.18p . for timelike (x' — x), if (x'° — x°) > the T-ordering will 
coincide with the P-ordering, while if (x'° — x°) < the T-ordering will be the opposite 
of the P-ordering (with a change of sign for the fermion fields), the color indices being 
already fixed from the P-ordering. We are in a situation which is very similar to the case 
of the ordinary two-point function, with the difference that for an n th -order term there 
are (n + 1) fields instead of two ((n — 1) gluon, one quark and one antiquark fields). 

Using for each of the two products which make the T-product the spectral analysis 
with intermediate states, taking into account the bounds on the parameters of the P- 
ordering and using causality, one arrives at a generalized form of the Kallen-Lehmann 
representation for the Green's function S in momentum space, in which the cut starts on 
the real axis from the quark mass squared m? and extends to infinity [271 1281 1291 1301 |3T] . 
The generalization is due to the fact that each gluon field is integrated along the path and 
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this introduces, when using for the latter a dimensionless parameter A varying between 
and 1, a multiplicative factor (x' — x), which is converted in momentum space into a 
derivation operator; each such factor increases by one unit the power of the denominator of 
the dispersion integral. Finally, because of the fact that we are dealing here with a gauge 
invariant quantity, we expect not to encounter at the end spurious infrared divergences. 

To summarize the above results, we introduce the Fourier transform of the Green's 
function S(x,x'), for which we also take into account translation invariance: 



S(p) has the following representation in terms of real spectral functions p™ and p$ 
(n = l,...,oo): 



This is a conservative representation of the various contributions encountered above; sim- 
plifications or recombinations into more compact forms might still occur. It is evident 
that formally the powers of the denominators can be lowered by integration by parts; 
however, possible singularities of the spectral functions at threshold could prevent such 
an operation. 

We assume that the above representation, obtained from the domain of perturbation 
theory, remains also valid in non-perturbative regimes. One expects that the resulting 
singularities are strong enough to screen the quark pole and other physical type singular- 
ities. 

Further study is needed to define more accurately the properties of the spectral func- 
tions. Nevertheless, representation (16. 21) . or another one related to it, might be tried for 
the investigation of the solutions of the corresponding integral equation. 

7 Summary and comments 

We have expressed the equation of motion of the gauge invariant quark two-point function 
having a straight line path as an integral or integro-differential equation involving the 
series of all two-point functions with paths of skew-polygonal type, in which the kernels 
are given by quark Green's functions and rigid path derivatives of the logarithms of the 
Wilson loop averages with contours made of these lines. 

Gauge invariant quark Green's functions also satisfy, in addition to their equations of 
motion related to the quark ends, other equations of motion resulting from local defor- 
mations of their paths. The latter equations are typically those of the path-ordered phase 
factors and once Wilson loops are introduced through the calculations, they reduce to 
the characteristic equations of Wilson loops, i.e., to the Bianchi identity and to the loop 






(6.2) 
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equation or Makeenko-Migdal equation [6j [7J [8] . We have not insisted on that aspect of 
the problem, since it has been widely studied in the literature. This implies that when 
the Wilson loop averages are used in the kernels of the above integral equations, they are 
understood as being solutions (at least approximately) of their own equations of motion. 

We have emphasized the fact that the series of kernels appearing in the integral equa- 
tion can be considered, on quantitative grounds, as a perturbation series simultaneously 
for short- and large-distances and therefore could be approximated, for a starting calcu- 
lation, by its lowest-order non- vanishing term. 

A question which was not considered in the present work concerns the renormalization 
properties of the gauge invariant Green's functions. These seem to be intimately related, 
through the integral equation, to those of the Wilson loop averages [HJ [15] and could be 
dealt with only when explicit expressions of the latter are introduced. 

The method of functional relations between two-point Green's functions with different 
numbers of segments along their paths, can also be applied, with appropriate general- 
izations, to other n-point Green's functions, with n > 4. That problem is particularly 
relevant for the derivation of a bound state equation for a quark- ant iquark system. 
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A Summation method 

Integral equations of ordinary Green's functions represent in general the result of summing 
the classes of reducible diagrams in terms of free propagators. One should expect that a 
similar procedure might also be operative in the case of gauge invariant Green's functions. 

To implement this method of approach, we should start with expressions involving 
free quark propagators. To this end, we have to use either of the two representations 
(13.41) or (13.51) of the quark propagator in external field and expand S(A) in the defining 
equation of 5* [Eq. (I4.2p ]. Each of the two representations has its own advantages and 
could be preferred for a definite aim. Thus, representation (13.41) is more appropriate to 
obtain rapidly the structure resulting from equations of motion with respect to x, while 
representation (I3.5P is more appropriate for the calculation of SS(x, x')/5x~. 

Using first representation f)3.4j) . the expansion of S(A) in Eq. (14. 2p generates Wilson 
loops with skew-polygonal type contours, accompanied with free quark propagators: 

S(x, x>) = S (x, x>) e F *&> x ) - S (x, yi ) 7 ai S ( Vl , x>) e F ^ x '> x ' 

oo 

+ S {x, yi ) 7 Q1 S (yi, m) l a2 ■■■ l a > S ( yj , x') 

J =2 
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A similar expansion can also be done for S(n) (n > 1), starting from Eq. 



S {n) (x, x'- z n . u ..., Zl ) = S (x, x') e F ^> z n-h • • • , *i, z) 
-S (x, Vl ) 7 ai S { yi , x') r 4 TT e F -+2(x', zn-u ...,z 1 ,x,y 1 ) 

oyi 

oo 

+ ^oO, yi) 7 Ql Sofoi, y 2 ) 7 ° 2 ■ ■ ■ 7° J x') 

J =2 

x j S ^_ p Fn+j+i( x ^ z n-u---,z 1 ,x,y 1 ,...,y j ) (A 9] 

X 5y? 1+ 5y? + 6y? + ' { ] 

Use of representation (13.51) yields equivalent expressions for S and S( n ) (n> 1): 

oo 

S(x, x') = £ S (x, Vl ) 7 ai S (yu y 2 ) 7° 2 ■ ■ ■ 7°' S {y j , x') 

3=0 

x = 5 - = 5 - _ 5 - e F j+2( x 'i x,yt,..., y 3 ) ? ( A3 ) 



S( n )(x, x'; z n _i, . . . , z x ) = S (x, yi) 7 Ql S (yi, y 2 ) l a2 ■■• l aj So(y 3 ,x') 

3=0 

x _J? _ !? ^_F n+i+1 (j/, a:, yi,...,yj) /a 4 n 



Equations (1A.3I) and (1A.4I) could also have been obtained from Eqs. (lA.lj) and (1A.2I) 



respectively, by integrations by parts; at the internal junction points 7/j of the segments 
of the paths, one has the equivalence relations d/dy i = 8/8yf + 6/5y^ [Eq. (12.101) ]. 

The equations of motion with respect to x can be evaluated easily from representations 
(lA.ip and (IA. 2|) . Because of the appearance of delta-functions from the propagators 
So(x,yi), there are cancellations between successive terms and one finds: 

(i~f.d(x) — m)S(x, x') = iS 4 (x — x') 

oo 

+ il" S (x, Vl ) 7 ai S ( yi , y 2 ) 7 Q2 ■ • • 7 Qj Sofa, x') 

3=0 

X J. I I d _ e F j+ 2(x',x, yi ,..., yj ) 

5x»- 5y^ + Sy? + 5y? + 

= i8\x-x') + iY ES j^'\ (A.5) 
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(i7-<9(*) - m)S (n) (x, x'; z n . u . . . , z x ) = i5\x - x') e F "( x ' Zn ~^ ■ ■ ■ ' z *> 

oo 

+ il" So(x, yi) 7 Q1 S {yi, y 2 ) 7° 2 ■ ■ ■ l a > x') 

3=0 



5 5 5 5 



?Fn+j+l( x ' i z n-\i ■ ■ ■ j %li x i Vli ■ ■ ■ j Vj) 



5x»-5y^5y? + 5y] 
= i5\x - x') z -i> ■ ■ ■ > *0 + ^W^^-i'-- 2 ') , (A . 6) 

Considering representation ( 1A.3I) . one immediately checks that it has the structure of 
the integral equation (15.91) ; this represents of course the integrated form of the equation 
of motion (1A.5I) . 

The action of the rigid path derivatives on the exponential functionals in Eqs. flA.ll) - 
(1A.4I) can be evaluated easily. The aim is then to group the various terms that appear in 
the expression of 5S(x, x')/5x~ to bring the latter into the form of Eq. ( 15.81) . We shall do 
this in a perturbative expansion with respect to the number of derivations, by retaining 
up to derivatives of third order and showing that Eq. (15. 8p can be obtained up to the 
terms. That approximation is sufficient to illustrate the various aspects of the method 
under consideration. 

For the calculation of 5S(x, x')/5x~ , it is preferable to start with representations ( 1A.3I) 



and (IA.4I) of S and S( n y, we shall indicate below the specific differences one meets when 
starting with representations (lA.lj) and (1A.2j) . although the final result is the same. The 



advantage of the former representations is that when a derivation 5/Sy^ acts on Su+i) in 
which yi is a fermionic end, it can be directly replaced in terms of the equation of motion 
operator and a delta-function, allowing an integration by parts. This is not the case with 
the operator 5/5yf, which first should be transformed into 5/5y^ before an integration 
by parts be possible. In the expansions (1A.3[) and (1A.4|) . S will be approximated with 



the first three terms, S(2) with the first two and S^) with the first term. We shall also 
assume the backtracking property of the path-ordered phase factors [10], which means 
that U(y,x)U(x,y) = 1, the same path being run forth and back. This means that 
W 2 = 1 and F 2 = 0. 

Calculating 5S(x,x')/5x~ from Eq. (1A.3|) one obtains: 



5S ( X ' X> ) = S r x x >) 6 e F 2 (x',x) , s r v A~aisJ v , x ') — —eM^^^yi) 



JLJ!_JL_ x, y u y 2 ) 



+S (x, yi ) 7 Q1 S (yu y 2 ) l a2 S (y 2 , x') ^_ „ 
+ ■■•• 1 2 (A.7) 

Calculating the derivative 5/5x~ , bringing the result to the left and completing in the 
first, second and third terms of the right-hand side of Eq. ( 1A.7I) the functions S, S(2) and 
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5(3), respectively, we obtain: 



6S(x,x') _ 5F 2 (x',x) 
8x^ Sx 11 - 
, 5 2 F 3 (x',x,yi; 



5(x,x') + s (x,m)7 Q1 



<5F 3 (o/,a:,2/i) 5F 2 (x',x) 



Sxv- 



S(2)(yi,x';x) 



+ ( 



5xf-8y^- 

S (x, yi h ai S (yi,y2h a2 
5 2 F, 



5F 4 (x',x,yi,y 2 ) 6F 3 (x',x,y 1 )\ 5 5 



5x^~ 



6 2 F d 



^Sx^-Sy" 1 Sx^Sy? 



5 2 F 4 



Sx^-Sy^ 2 ' 



5yr 



+ 



5 3 F 



5x»-5 y ^-5y% 2 - 



5x^~ 



S(3)(y 2 ,x';x,y 1 ). 



ty? 1 ' Syl 



a 2- 



(A. 



(Higher-order terms in the derivatives are neglected.) The operators S/Sy" 1 * and S/5y 2 2 ~ 
acting on 5(2) and 5(3), respectively, can be replaced in terms of the corresponding equation 
of motion operators and delta functions and then integrations by parts are carried out; 
simailarly, the operator S/Sy" 1 ' acting on 5(3) can be replaced by d/dy 1 — S/Syf followed 
by an integration by parts. One finds at the end: 



6S(x,x') _ SF 2 (x',x) 
5x^~ 8x^~ 
5 3 F 4 (x',x,y 1 ,y 2 ) 

5x»-8y? 1+ ty% 
6 2 F 4 5F 4 



S(x, x') 



6 2 F 3 (x',x, yi ) 



S (x,yi)>y ai S {2 )(yi,x';x) 



+ 



Sx^'Syi 1 ' 

h S Q (x, y x ) 7 ai S ( yil y 2 )~f a2 S {3) (y 2 ,x';x,yi) 



5x»-5y¥ + Sy^ 1 



:5 (x,?/i)7 Q!1 S (yi,y 2 )>y a2 S {3) (y 2 ,x'; x,y x ). 



(A.9) 



This result could also have been obtained by using the method of Sees. HlandEl but using 
for the expansion of S(A) Eq. (13.51) . instead of (I3.12p . It is sufficient for this to replace 
in Eq. (15.81) 55/5^ by zero and 5 in internal lines by So. 

If we had started for the previous calculations from representations ( 1A.1I) and ( 1A.2I) . 
we would have found in a first stage the terms of Eq. ( 1A.9I) with a remainder con- 
taining difference terms of the type (5 2 F 4 /5x~5yi — S 2 F 3 /Sx~Syi)S/5y 2 e Fi , (5F 4 /5x~ — 
5F 3 /5x~)S 2 /SyfEy^e^, etc. Integrations by parts show that the remainder is null. For 
this, one must proceed in two steps. First one converts S/5y 2 into S/5y 2 by the formula 
S/Sy£ = 9/dy 2 — 5/5y 2 and an integration by parts is done with respect to the total 
derivative of y 2 . Second, for the term containing 5/5y 2 , one completes the exponential 
function with the multiplicative propagator So(y 2 ,x') into 5(3) and then the operator 
S/5y 2 is replaced in terms of the equation of motion operator and a delta function and 
a new integration by parts is done. The net result is zero. This method of calculation is 
repeated for all parts of the remainder. One thus finds the same result (1A.9I) from both 
representations (lA~ll - (lA~2|) and (lA~3l) - (IA~4l) . 
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Inspection of Eq. (1A.9I) shows that the last term is of the nested type, with 5F^/5yi 
representing a kind of self-energy insertion on the line xy 2 . It should naturally be grouped 
with the free propagator So appearing in front of S^) (with a relabelling of the variables 
y 1 and y 2 ) to produce the full Green's function S [Eq. ( 1A.1I) ] (at the present level of 
approximation). Nevertheless, we are faced with the phenomenon of the Wilson loop 
background effect: the various factors that appear in the last term of Eq. (1A.9I) involve 
_F 4 and not F 3 which is the required function in the next-to-leading term of S. The 
recombinations that we can do to reconstruct full Green's functions in internal lines leave 
at the end remainders of the nested type. 

Before proceeding to a recombination of the above factors in the general case, let us 
consider first the particular case of two-dimensional QCD in the quenched approximation 
(quark loops neglected) [13]. We assume that the Wilson loop contours that mainly 
contribute to the internal integrations are simple and convex, in particular without self- 
intersections. In that case, the logarithm of the Wilson loop average is given by the 
area of the surface delimited by the closed contour and the functions F$ are proportional 
to such areas. Furthermore, it is evident that these areas are separable into smaller 
ones. For the specific case above, we have F^(x',x,yi,y 2 ) = F 3 (x',x,y 2 ) + F 3 (x, yi, y 2 ), 
which also implies 5F A (x',x,y 1 ,y 2 )/Syf = 5F 3 (x,y 1 ,y 2 )/Syf, 5 2 F A (x\ x, y u y 2 )/5x~5y£ = 
5 2 F 3 (x',x,y 2 )/Sx~5y 2 , etc. With such decompositions, one easily transforms the last 
term of Eq. flA.9j) into a form that is absorbed by the term containing Sf 2 ) to yield in 
front of it the full Green's function S. Replacing in the remaining term containing S/g) 
the free propagators So by S (valid at the present level of approximation), one finds 

SS(x,x') 5F 2 (x',x) , 5 2 F 3 (x',x,y 1 ) , 

= S(x, x ) - S(x, yi ) 7 S {2) ( yi , x ; x) 

+ s Z% x ^tr S ^ Vl) 7 " s{m ' m) 7 " 5(3)Q/2 ' x>] x ' (A,10) 

which is an expansion with full Green's functions in internal lines and kernels of the 
irreducible type without nested diagrams. 

In four dimensions, the above decompositions of the Wilson loop averages are not 
generally valid and one has to evaluate the remainder with respect to Eq. (1A.10I) . To 
this end, we complete in Eq. (1A.9I) the factor S (x, y\) in the term containing S( 2 ) into 
S(x,yi) and isolate the rest, which is equal to 

- z.h^otz,, V, f 3 f£; Zl) eK<.yu*, *.) 7 - ] S m ( Vl , X <; ,) (A.ll) 

and which we write in the form 

-S (x, z 1 )^S (z 1 ,y 1 ) 5F *f£ Zl) e F ^ x > - y 2 ) 



5 2 F 4 {x\x,yi,y 2 ) . , . 
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Writing the delta-function in the form iS 4 (yi — y 2 ) = (ij-d( yi ) — m)So(yi, y 2 ) and making 
an integration by parts with respect to y\ and neglecting higher-order derivative terms, 
we obtain 

Q( \ aic/ jF 2 (yi,x) F 2 (y u x) a 2 ^ F ^W, x , Vu V2) Q , , , r\io\ 
-So(x,y 1 )j ai So(yi,y2) ^ ^ e 2W ' h — Sx»-6y% 2+ — (3) ^ 2 ' x ' x ' ^ ^ 

In the first terms we recognize the dominant pieces of the product ^y£p^- j ai S(yi,y 2 ). 
Using the equation of motion of S and making again an integration by parts with respect 
to y 1 (neglecting higher-order derivatives) we find the final expression 

cr ^ aJ S {VuV2) a J 2 F 4 (x',x,y u y 2 ) , , . r\-i A \ 
S(x, Vl h ^ ai - 7 ^-^+ g (3)(y 2 ^;x, yi ), (A.14) 



which, when grouped with the last term of Eq. flA.9j) . in which the Sqs may be replaced 



by Ss, yields the nested piece of Eq. 05. 8 p accompanying S^y. 

S 2 F A (x',x,y 1 ,y 2 ) ( 5S(y ll y 2 ) . , 5F 4 \ , 

5 x n-S y ?+ S &y^ ( ^1- + g (^i^2)^ IT j7 2 5 ( 3)(y 2 ,x;x, yi )- (A.15) 

We thus recover, together with the terms of Eq. flA.101) . the first terms of Eq. (15.81) . up 
to S*(3). The calculation could be continued to higher-orders in the derivative terms, but 
it rapidly becomes complicated and loses interest for practical applicability. The method 
is useful for low-order perturbative calculations and for analysis of general qualitative 
properties. 

B Quark self-energy 

In order to construct the quark self-energy function directly, without having recourse to 
the explicit expression of the Green's function S, we start from its relationship with the 
derivative terms of S. The self-energy function, which we designate by S, is defined from 
the inverse of the Green's function: 

iS-\x,x') = [i 1 .d {x) -m)5\x -x') - E(x,x'). (B.l) 

Comparison with the equations of motion (I3.8P and (13.91) yields: 

£(*, x") S(x", x') = vf S(x, x") E(x", x') = -i 6 -^^ Y- (B.2) 

Letting S 1-1 act on these equations, one ends up with the equation for S: 

£(*, x') = -vf 0^ l v + tT(x, y) S(y, y') E(y', x'). (B.3) 
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+ig 



The second-order derivative S 2 S(x, x')/5x ll ~5x u+ corresponds to a generalization of the 
first-order derivatives 5S/5x ± , defined in Eqs. (12.201) and (I2.6p - (I2.7I) . where now the two 
derivations act on the same segment xx'; in this case, one has also to take into account 
the contributions coming from coincident points. Explicitly, one has: 

f^fj + = (fc) a [ d\d\'(l - A) A' [17(1, A')^(A / )F /3l ,(^(A'))t/(A', A) 
xz' a (X)F a ,(z(X))U(X, 0) + U(l, \)z' a (\)F a ,(z^))U(\ A , )^(A , )F /3l/ ( 2 (A , ))f/(A', 0)" 
f d\{\ - A)A U(l, X)z' a (X) (V u F^(z(X)))u(X, 0) 

+ig f 1 d\(l - A) U{1, X)F^(z(X))U(X, 0), (B.4) 
J o 

where z'(X) = || = x' — x and V is the covariant derivative, (VF) = (dF) + ig[A,F]. 
The calculation was done by first deriving with respect to x and then with respect to x' . 
Had we interchanged the orders of derivation, the last two terms would be modified in 
the following way: z' a V ' w F afl — > z' a V ^F^, (1 — X)F UfM — > —XF^. The two expressions 
are, however, equivalent, since the orders of derivation are irrelevant, due to the Bianchi 
identity satisfied by F. This can be checked directly by taking the difference of the two 
expressions. 

The two derivative term ( IB. 41) requires a careful treatment, since it contains diver- 
gences or singularities not present in one derivative terms. This is the case for the trace 
of the tensor flB.41) : for neighboring points in the expression inside the brackets, the two 
Fs lead to a divergence, even when short-distance perturbative interactions are ignored 
[321 [26] : still for the trace, the term with the covariant derivative reduces to the gluon 
equation of motion operator and hence yields a delta- function (plus a quark current); 
these singular terms should be grouped together to set up a regularized form of the corre- 
sponding quantities. Concerning the traceless part of the coincident points contribution 
(the last two terms of Eq. ( 1B.4I) ). we observe that it is of order g and not g 2 ; this implies 
that when expanding equation ( IB. 31) in terms of derivatives of Wilson loop averages, one 
should count the latter contribution as a one derivative term. Finally, the role of the last 
term of Eq. (1B.3[) is to cancel similar reducible type terms that might emerge from the 



expansion of the second-order derivative piece. Nevertheless, because of existing back- 
ground effects, as in the case of nested diagrams met in Sec. El the cancellations are only 
partial. 

The above features make the direct treatment of the self-energy function rather intri- 
cate and less appealing than that of the Green's function itself. This underlines the fact 
that proper vertices do not seem to play a primary role in the present approach. We shall 
not pursue any longer here the study of the self-energy function. 
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